PB96-148804 


June  1993 


Report  No.  STAN-CS-93-1478 


PB96-148804 


Two-Handed  Assembly  Sequencing 


by 


R.H.WiIson,  L.  Kavraki,  T.  Lozano*Perez,  J.C.Latombe 


Department  of  Computer  Science 

Stanford  University 
Stanford,  California  94305 


REPRODUCED  BY:  NTtS 

U.S,  Department  of  Commerce 
National  Technical  Information  Service 
Springfield,  Virginia  22161 


BIBLIOGRAPHIC  INFORMATION 


PB96-148804 


Report  Nos:  STAN-CS-93-1478 

Title:  Two-Handed  Assembly  Sequencing. 

Date:  cJun  93 


Authors :  R.  H.  Wilson,  L.  Kavraki ,  T.  Lozano-Perez,  and  J.  C.  Latombe. 

Performing  Organization:  Stanford  Univ^,  CA.  Dept,  of  Computer  Science. **Massachusetts 
Inst,  ot  Tech!,  Cambridge.  Artificial  Intelligence  Lab. 

Sponsoring  Organization:  ^Defense  Advanced  Research  Projects  Agency,  Arlington, 
\/A.*Uffice  OT  Naval  Kesearch,  Arlington,  VA. 

Contract  Nos :  ONR-N00014-92-J-1809,  ONR-N00014-91-J-4038 

NTIS  Field/Group  Codes:  41C  (Robotics/Robots) ,_41E  (Manufacturing,  Planning, 

Processing  &  Lontrol JT  62  (Computers,  Control  &  Information  Theory),  94G 
(Manufacturing  Processes  &  Materials  Handling) 

Price:  PC  A03/MF  AOl 

Availability:  Available  from  the  National  Technical  Information  Service,  Springfield, 
VA.  22i6i 

Number  of  Pages:  29p 

Keywords:  *Robotics.  *Algorithms.  *Assembling,  Geometry.  Components.  Mathematics, 
Automat 1  on ,  Matching.  Fitting,  Computer  programs.  Computations,  Assembly  processes. 

Abstract:  This  paper  considers  the  computational  complexity  of  automatically 
determining  assembly  sequences  for  mechanical  products.  Specifically,  we  address  the 
partitioning  problem:  given  an  assembly  of  rigid  parts,  identify  a  proper  subassembly 
that  can  be  removed  as  a  rigid  object  without  disturbing  the  rest  of  the  assembly.  We 
examine  the  complexity  of  the  partitioning  problem  under  varying  types  of  relative 
motions  allowed  for  the  subassemblies.  We  show  that  when  arbitrary  motions  are  allowed 
to  separate  the  two  subassemblies,  partitioning  is  NP-complete.  We  then  describe  a 
general  framework  for  reasoning  about  assembly  motions  called  the  interference 
diagram  In  its  most  general  form  the  interference  diagram  yields  an  exponential -time 
algorithm  to  partition  an  assembly.  However,  two  special  cases  of  the  interference 
diagram  studied  in  this  paper  yield  polynomial -time  sequencing  algorithms  The  first 
case  occurs  when  assembly  motions  are  restricted  to  single  translations.  The  second 
case  considers  infinitesimal  rigid  motions  in  translation  and  rotation,  and  yields  a 
superset  of  all  feasible  partitionings. 
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1  Introduction 


Automated  assembly  of  composite  products  has  been  a  goal  of  robotics  researchers 
since  the  beginnings  of  the  field  [9,  17],  and  robots  are  used  today  to  assemble  a 
variety  of  products.  However,  programming  these  robots  remains  tedious  and  error- 
prone,  and  algorithms  are  sought  to  help  program  them.  Research  in  assembly 
planning  aims  to  automate  or  partially  automate  the  ordering  of  assembly  operations 
and  selection  of  the  required  tools  and  motions,  as  well  as  provide  fast  manufacturing 
evaluation  of  assembly  designs.  This  paper  concerns  a  subproblem  of  assembly 
planning  called  assembly  sequence  planning,  or  assembly  sequencing,  in  which  the 
relative  motions  of  the  parts  of  the  assembly  are  determined  without  considering 
the  tools,  fixtures,  or  robots  required  to  achieve  those  motions. 

Disassembly  planning  has  been  a  very  popular  approach  to  assembly  sequenc¬ 
ing  (see  for  instance  [1,  4,  7,  15,  16,  19,  20,  28,  31]).  The  most  constraints  on 
assembly  motions  are  present  in  the  final  assembled  state  of  the  product,  so  reason¬ 
ing  backward  from  the  assembled  state  reduces  the  branching  factor  of  the  search 
considerably.  The  result  is  a  sequence  of  relative  motions  for  the  parts  called  a  dis¬ 
assembly  sequence.  Here  we  assume  the  parts  are  rigid,  so  a  disassembly  sequence 
is  the  reverse  of  an  assembly  sequence,  and  conversely. 

In  addition,  we  limit  consideration  to  monotone  binary  assembly  sequences: 
those  that  at  each  step  merge  exactly  two  rigid  subassemblies  to  make  a  larger 
one.  Such  sequences  are  also  called  monotone  two-handed,  since  one  hand  is  needed 
to  hold  each  rigidly  moving  subassembly  (a  table  or  fixture  counts  as  a  hand,  if 
one  is  used)  [23].  For  example,  the  assembly  in  figure  1(a)  can  be  assembled  by 
a  monotone  binary  assembly  sequence,  while  the  assembly  in  figure  1(b)  cannot; 
the  latter  requires  one  part  to  be  placed  in  an  intermediate  position  then  moved 
again  later.  This  limitation  is  well  supported  by  typical  assemblies  seen  in  industry, 
and  most  experimental  assembly  planning  systems  today  find  only  monotone  binary 
sequences. 

The  success  of  such  a  disassembly  planner  depends  on  solving  the  partitioning 
problem:  given  an  assembly,  identify  a  proper  subassembly  that  can  be  removed  as  a 
rigid  object  without  moving  the  rest  of  the  parts.  The  two  subassemblies  produced 
can  be  partitioned  in  turn,  and  so  on  until  only  individual  parts  remain.  As  a  result, 
the  partitioning  problem  is  the  core  of  the  monotone  binary  assembly  sequencing 
problem. 

In  this  paper  we  consider  the  computational  complexity  of  the  partitioning  prob¬ 
lem.  We  show  that  when  the  removed  subassembly  may  foUow  an  arbitrary  path 
(or  an  arbitrary  sequence  of  translations),  partitioning  is  NP-complete,  even  for 
two-dimensional  assemblies  of  polygons.  We  then  introduce  a  general  framework  for 
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Figure  1:  (a)  is  a  monotone  binary  assembly  while  (b)  is  not. 


reasoning  about  the  motions  of  subassemblies  called  the  interference  diagram.  The 
interference  diagram  makes  expUcit  the  constraints  on  motion  of  all  parts  of  the 
assembly  relative  to  all  others,  and  embodies  the  monotone  and  binary  restrictions 
exactly.  Under  this  framework  we  analyze  the  conditions  under  which  partitioning 
is  tractable.  While  the  interference  diagram  for  an  assembly  can  be  constructed  in 
polynomial  time,  in  its  most  general  form  it  yields  a  partitioning  algorithm  that  is 
exponential  in  the  number  of  surfaces  describing  the  assembly. 

However,  restrictions  on  the  motions  allowed  to  separate  two  subassemblies  give 
rise  to  special  cases  of  the  interference  diagram  that  can  be  analyzed  in  polynomial 
time.  We  analyze  two  cases  here.  In  the  first  case,  assembly  motions  are  restricted 
to  single  translations  to  infinity,  and  the  resulting  algorithm  identifies  a  removable 
subassembly  in  time,  for  an  assembly  of  n  polyhedra  having  a  total  of 

V  vertices.  In  the  second  case,  infinitesimal  rigid  motions  (including  translation 
and  rotation)  are  considered;  the  set  of  movable  subassemblies  thus  identified  is 
a  superset  of  those  that  can  be  removed  by  arbitrary  paths.  A  subassembly  that 
can  translate  or  rotate  a  small  distance  from  the  rest  of  the  assembly  is  identified 
in  0(c®  +  m^c^)  time,  for  a  polyhedral  assembly  with  m  pairs  of  parts  in  contact, 
whose  contacts  are  described  by  c  €  O(v^)  point-plane  contact  constraints. 

2  Related  Work 

While  we  consider  only  monotone  binary  assembly  sequences  in  this  paper,  other  au¬ 
thors  have  examined  more  general  cases.  In  non-monotone  assembly  sequences,  parts 
may  assume  intermediate  positions  in  subassemblies.  Natarajan  [23]  and  Wolter  [31] 
show  that  assembly  sequencing  is  PSPACE-hard  when  non-monotone  sequences  are 
allowed.  Hoffman  [14]  describes  one  of  the  few  non-monotone  assembly  planning 
systems. 

In  non-binary  assembly  sequences,  more  than  one  subassembly  may  move  in- 
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dependently  at  the  same  time.  Palmer  [25]  considers  the  infinitesimal,  non-binary 
partitioning  problem:  determining  whether  a  feasible  set  of  simultaneous,  infinitesi¬ 
mal  motions  exists  for  the  parts  of  a  polygonal  assembly.  He  shows  that  this  problem 
is  NP-complete  by  a  reduction  from  3-Satisfiability.  However,  few  industrial  assem¬ 
blies  require  the  simultaneous  independent  motion  of  many  parts.  In  this  paper  a 
subassembly  always  refers  to  a  set  of  parts  in  fixed  relative  position. 

Most  experimental  assembly  planning  systems  find  only  monotone  binary  as¬ 
sembly  sequences.  These  systems  address  the  partitioning  problem  in  a  variety  of 
ways.  Several  systems  [4,  16,  20]  enumerate  all  possible  subassemblies  and  test  each 
for  removal,  an  approach  that  is  practical  for  small  assemblies  but  exponential-time 
in  the  worst  case.  Another  approach  is  to  consider  only  sequences  in  which  each 
operation  mates  a  single  part  with  a  subassembly  [19,  30].  This  approach  reduces 
the  partitioning  problem  to  checking  for  the  motion  of  each  part,  but  risks  losing 
some,  maybe  all,  of  the  assembly  sequences  for  a  product.  Finally,  many  systems 
predefine  a  small  set  of  assembly  trajectories,  or  heuristically  choose  trajectories 
based  on  features  of  the  parts  or  their  contacts  [3,  7,  14,  31].  However,  this  ap¬ 
proach  requires  extra  planning  information  be  present  in  the  assembly  description, 
and  risks  missing  some  assembly  sequences. 

Arkin,  Connelly,  and  Mitchell  [1]  use  the  concept  of  a  monotone  path  among 
obstacles  to  derive  a  polynomial-time  algorithm  for  partitioning  an  assembly  of 
polygons  in  the  plane  with  a  single  infinite  translation.  This  corresponds  to  the  first 
special  case  of  the  interference  diagram  considered  below.  However,  their  approach 
is  significantly  different  from  ours  and  does  not  directly  extend  to  the  3D  case. 

In  the  rest  of  this  paper,  we  analyze  the  partitioning  problem  from  a  purely  geo¬ 
metric  and  computational  point  of  view.  In  section  3  we  show  that  when  an  arbitrary 
translational  path  is  allowed  to  separate  two  subassembbes,  the  partitioning  problem 
for  assemblies  of  simple  polygons  in  the  plane  is  NP-complete.  Section  4  introduces 
the  interference  diagram,  a  novel  configuration-space  formulation  of  the  partition¬ 
ing  problem.  W  hile  the  interference  diagram  for  an  assembly  can  be  constructed  in 
polynomial  time,  its  analysis  requires  exponential  time.  However,  in  sections  5  and  6 
we  show  that  when  assembly  motions  are  restricted  to  single  translations  and  in¬ 
finitesimal  rigid  motions,  respectively,  the  resulting  special  cases  of  the  interference 
diagram  yield  polynomial-time  partitioning  algorithms.  Finally,  section  7  examines 
the  boundary  between  tractable  and  intractable  cases  of  the  partitioning  problem. 


3  Complexity  of  Partitioning 

The  assembly  partitioning  problem  consists  of  identifying  a  subassembly  of  a  given 
assembly  that  can  be  removed  (as  a  rigid  object)  without  disturbing  the  other  parts 
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of  the  assembly.  In  this  section  we  show  that  deciding  whether  such  a  subassembly 
exists  is  an  NP-complete  problem.  We  outline  the  proof  for  an  assembly  of  polygons 
that  can  translate  and  rotate  in  the  plane.  A  more  detailed  version  of  the  proof 
can  be  found  in  [18].  The  NP-completeness  result  also  extends  to  some  restricted 
variants  of  the  partitioning  problem  that  are  presented  at  the  end  of  this  section. 

Let  a  rigid  motion  of  a  subassembly  5  be  a  set  of  simultaneous  motions  (trans¬ 
lations  and  rotations)  of  the  parts  of  S  that  preserve  the  relative  positions  of  these 
parts  throughout  the  motion.  The  subassembly  S  is  then  called  a  rigid  subassembly. 

Planar  Partitioning  (PP)  Given  a  set  A  of  non-overlapping  polygons  in  the 
plane,  decide  if  there  is  proper  subset  S  of  A  that  can  be  separated  from  A\S  by  a 
collision-free  rigid  motion  of  S. 

We  will  show  that  PP  is  NP-complete.  It  is  clearly  in  NP,  since  a  nondetermin- 
istic  algorithm  can  guess  S  and  then  invoke  a  path  planner  to  find  the  path  of  S 
out  of  the  assembly.  Schwartz  and  Sharir  [27]  have  shown  that  path  planning  can 
be  done  in  polynomial  time  for  the  case  considered  here. 

We  show  that  PP  is  NP-hard  by  a  reduction  from  .3- Satisfiability  (3-SAT),  a 
well  known  NP-complete  problem  [10].  An  instance  of  3-SAT  is  a  set  of  clauses 
C  =  {ci,  C2, . . . ,  Cm}  on  a  set  of  boolean  variables  U  =  {ui,U2, . .  - ,  u„},  where  each 
clause  is  a  disjunction  of  3  terms.  A  term  is  either  a  variable  u;  or  the  negation  of 
a  variable  ui.  The  problem  is  to  determine  if  there  exists  a  truth  assignment  of  the 
variables  that  satisfies  the  conjunction  of  the  clauses. 

For  any  instance  of  3-SAT,  we  construct  in  polynomial  time  an  assembly  of  non¬ 
overlapping  polygons  that  can  be  partitioned  iff  a  satisfying  truth  assignment  exists. 
Figure  2  shows  the  outside  box  and  the  key  of  the  constructed  assembly.  Other  parts 
are  contained  in  the  assignment  mechanism  and  the  AND/OR  mechanism,  detailed 
in  figures  3  and  4  respectively.  Our  construction  is  summarized  in  the  following: 

•  A  rigid  subassembly  must  be  removed  to  partition  the  assembly  A.  Call  it  5 . 

•  The  part  labeled  as  key  in  figure  2  blocks  the  only  exit  gate  of  the  assembly. 
Hence,  it  must  be  a  member  of  S. 

•  The  key  can  be.  removed  only  through  the  assignment  construct.  For  this  to 
happen  some  other  parts  of  the  assignment  construct  must  be  in  S.  These 
parts  represent  a  truth  assignment  for  the  variables  of  the  3-SAT  instance  arid 
move  rigidly  with  the  key. 

•  5  can  be  removed  only  through  the  AND/OR  mechanism.  This  mechanism 
enforces  the  clauses  of  the  3-SAT  instance. 
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Figure  2:  A  sketch  of  the  final  assembly. 


Figure  3:  The  assignment  construct  (not  drawn  to  scale). 


*  Finally  S  is  moved  to  the  cavity  on  the  upper  left,  rotated,  and  removed 
through  the  gap  left  by  the  key. 

The  assignment  mechanism  (figure  3)  consists  of  (i)  the  walls  of  the  assignment 
that  are  drawn  in  grey,  (ii)  the  key  which  is  a  3  x  1  rectangle  drawn  in  black,  and 
(iii)  two  3x1  white  rectangles  for  each  of  the  variables  of  the  3-SAT  instance.  One 
of  the  two  rectangles  that  correspond  to  the  variable  Uj  is  labeled  with  £/),  and 
the  other  with  Uj.  Notice  that  Ui  is  placed  always  on  top  of  TTi  in  the  assignment 
construct.  If  Uj,  {Ui  resp.)  is  a  member  of  5,  we  consider  that  the  truth  assignment 
true  {false  resp.),  has  been  chosen  for  the  variable  Ui.  We  indicate  with  xq  the  initial 
position  of  the  key  and  with  a:,-,  the  initial  position  of  the  assignment  rectangles  for 
the  variable  Ui,  i  =  l,...,n.  The  choice  of  the  x,’s  is  crucial  and  it  is  described 
below. 

From  figure  2  it  is  clear  that  the  key  initially  can  move  only  to  the  right.  We 
observe  that  the  key  will  not  be  able  to  pass  through  the  assignment  mechanism  if 
no  other  parts  of  the  mechanism  are  moved.  In  addition,  the  parts  removed  must 
translate  rigidly.  Hence,  the  only  subassembly  5  that  stands  a  chance  to  move 
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out  of  the  assignment  mechanism,  and  eventually  out  of  the  total  assembly,  is  a 
subassembly  that  consists  of  the  key  and  at  least  one  of  Ui  or  Z7,-,  for  each  i.  Since 
the  exit  gate  of  the  assignment  mechanism  has  a  height  of  only  1,  exactly  one  of 
Ui  or  Ui,  for  each  i,  must  be  selected  in  S.  Let  L  be  the  length  of  the  moving 
subassembly  5. 

The  coUision-free  motion  of  S  out  of  the  assignment  construct  is  possible  only  if 
the  Xi's  are  selected  carefully.  We  show  in  the  appendix  how  to  choose  the  Xi’s  to 
guarantee  the  following;  for  any  k,  when  Uk  or  Uk  (whichever  of  them  belongs  to 
S)  is  close  to  position  Xp,  p  k,  and  needs  to  go  through  the  hole  that  has  been 
created  at  this  position,  none  of  the  other  assignment  rectangles  or  the  key  is  close 
to  a  position  Xi,  i  7^  p.  In  other  words,  aU  the  other  parts  of  S  are  in  the  wide 
free  sections  of  the  assignment  mechanism  and  can  follow  the  constrained  motion 
of  the  part  that  is  close  to  Xp.  Hence,  the  above  mechanism  ensures  independent 
selection  of  the  variable  assignments  of  the  3-SAT  instance  and  allows  the  resulting 
subassembly  to  translate  out  of  the  construct. 

Once  S  translates  out  of  the  exit  gate  of  the  assignment  construct  it  must  pass 
through  the  AND/OR  mechanism.  This  mechanism  is  a  sequence  of  OR  gates,  one 
for  each  clause  of  the  3-SAT  instance.  The  OR  gate  for  the  clause  V  uj  V  Uk  is 
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shown  in  figure  4.  We  observe  from  the  figure  that  there  are  three  possible  ways  for 
S  to  go  from  GateA  to  GateB.  Each  of  these  enforces  a  truth  assignment  for  one 
of  the  terms  of  the  clause.  Suppose  for  example  that  5  goes  through  the  section 
enclosed  in  the  dashed  box  in  figure  4.  This  section  enforces  the  truth  assignment 
true  for  the  variable  n,-.  Here  is  why:  when  the  key  is  at  xo,  the  rectangle  chosen 
for  the  truth  assignment  of  u,-  is  at  position  x,-.  Unless  Ui  is  selected  in  S,  it  is 
impossible  for  S  to  go  through  the  dashed  box  of  figure  4.  Notice  however  that  the 
rest  of  5  can  be  threaded  through  the  gates  at  xo  a-nd  x,  without  problems:  because 
of  the  property  of  the  x,’s  mentioned  above,  when  a  part  of  S  needs  to  go  through 
the  above  gates,  none  of  the  other  parts  of  S  is  close  to  a  narrow  passage.  Hence,  5 
can  follow  the  motion  of  its  constrained  part  without  being  obstructed  by  the  walls 
of  the  OR  gate. 

Suppose  there  exists  a  satisfying  truth  assignment  for  the  3- SAT  instance.  Let 
5  be  the  subassembly  that  consists  of  the  key  and  encodes  this  truth  assignment. 
S  can  go  through  the  AND/OR  mechanism  since  it  can  translate  through  each  of 
its  OR  gates.  Then  S  can  translate  to  the  upper  left  corner  of  the  assembly,  rotate 
there  by  90  degrees  and  exit  through  the  2-unit  wide  gate  that  was  initially  blocked 
by  the  key. 

Conversely,  assume  that  the  assembly  in  figure  2  can  be  partitioned  and  let  S 
be  the  subassembly  that  is  removed  from  it.  5  clearly  contains  the  key.  As  we 
argue  above,  the  key  can  be  removed  only  in  a  subassembly  that  contains  a  truth 
assignment  for  the  variables  of  the  3-SAT  instance.  Since  S  can  pass  through  the 
AND/OR  mechanism,  it  represents  a  satisfying  truth  assignment  for  the  3-SAT 
instance.  Finally,  it  can  be  shown  easily  that  the  reduction  presented  above  is 
polynomial  in  the  size  of  the  3-SAT  problem.  It  follows  that: 

Theorem  1  PP  is  NP-complete. 

Variants  of  the  above  proof  establish  that  the  following  problems  are  NP- 
complete  [18,  29]: 

•  Partitioning  of  a  planar  assembly  when  rotation  is  not  allowed. 

•  Partitioning  of  an  assembly  of  polyhedra  in  space  when  rotation  is  allowed 
and  when  rotation  is  not  allowed. 

•  Partitioning  of  assemblies  with  parts  of  constant  complexity,  i.e.  each  part  is 
limited  to  a  constant  number  of  vertices. 

•  Partitioning  of  assemblies  of  constant  complexity  polyhedra  with  the  addi¬ 
tional  requirement  that  the  two  subassemblies  produced  are  connected.  A 
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subassembly  is  connected,  when  the  union  of  the  spaces  occupied  by  its  parts 
is  connected. 

•  Partitioning  of  assemblies  that  can  be  fully  disassembled.  The  reader  may  note 
that  the  assembly  in  figure  2  cannot  be.  This  result  proves  that  monotone  two- 
handed  assembly  sequencing  is  an  NP-complete  problem. 

4  The  Interference  Diagram 

In  this  section  we  present  a  general  framework  for  reasoning  about  assembly  motions, 
called  the  interference  diagram.  The  interference  diagram  is  motivated  by  the  desire 
to  reason  about  a  trajectory  for  a  subassembly  S,  without  knowing  a  prion  which 
parts  of  the  assembly  are  members  of  5.  Once  a  path  is  known,  or  as  it  is  constructed, 
constraints  on  which  parts  of  the  assembly  may  be  in  5  are  analyzed  to  determine 
the  members  of  S. 

The  interference  diagram  is  a  somewhat  complicated  configuration-space  dia¬ 
gram,  so  we  can  illustrate  it  graphically  only  in  simple  examples  such  as  the  one  in 
figure  5(a).  To  simplify  the  presentation,  we  will  limit  ourselves  to  pure  translations. 
Furthermore,  we  have  drawn  the  parts  with  substantial  clearances  so  that  each  of 
the  relevant  regions  in  configuration  space  would  be  full-dimensional.  We  sketch  the 
generalization  to  higher  dimensions,  rotations,  and  situations  involving  contacts  in 
section  4.4. 

4.1  Placing  a  Part 

Since  we  do  not  know  which  parts  wiU  be  moving  and  which  are  stationary,  the 
interference  diagram  represents  the  constraints  on  motion  between  any  pair  of  parts 
in  the  assembly.  For  every  pair  (X,  T)  of  parts,  the  set  of  placements  of  X  in  which 
it  intersects  with  Y  is  the  configuration  space  (or  C-space)  obstacle  for  the  pair  of 
parts  [21].  The  C-space  obstacle  for  part  X  moving  with  Y  as  an  obstacle  is  labeled 
X/Y,  while  the  C-space  obstacle  for  Y  moving  with  X  as  an  obstacle  is  labeled 
Y/X.  In  translation,  the  C-space  obstacle  X/Y  is  given  by 

X/Y  =  Y  e  X  =  {y  -  X  \  X  e  X,  y  e  Y} 

i.e.,  the  Minkowski  difference  of  the  two  sets  of  points  Y  and  X .  Note  that  Y/X  is 
simply  X/Y  rotated  by  x  radians  or,  equivalently,  every  forbidden  translation  x  of 
X  has  a  corresponding  forbidden  translation  —x  of  T. 

Figure  5  shows  the  6  C-space  obstacles  for  the  simple  example  we  are  considering. 
A  crucial  point  about  these  C-spaces  is  that  they  are  all  constructed  using  the  same 
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Figure  5:  C-space  diagrams  for  each  pair  of  parts  in  the  assembly.  The  notation 
X/y  indicates  that  X  (moving)  collides  with  Y  (stationary)  in  that  region. 
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reference  pointy  indicated  by  a  solid  circle  in  the  figure;  the  reference  point  of  an 
object  need  not  be  inside  that  object.  That  is,  treat  this  point  as  being  rigidly 
attached  to  each  part  in  turn.  Then  the  C-space  obstacles  represent  the  positions 
of  this  reference  point  for  which  two  parts  collide. 

To  remove  one  part  from  the  assembly,  these  pairwise  C^spaces  clearly  give  us 
the  information  we  need  to  plan  the  motion.  We  can  test  if  part  5,  for  example,  can 
be  moved  by  examining  whether  a  motion  exists  in  the  region  outside  the  union  of 
B/A  and  B/C.  However,  we  wish  to  reason  about  the  motion  of  subassemblies.  In 
principle,  we  could  enumerate  all  subassemblies  and  compute  the  C-space  obstacles 
for  them,  but  there  is  an  exponential  number  of  subassemblies  to  consider. 

4.2  Placing  a  Subassembly 

Instead,  we  reason  about  the  motion  of  an  unknown  rigid  subassembly,  while  gather¬ 
ing  constraints  on  which  parts  can  follow  that  motion.  Since  the  parts  of  the  moved 
subassembly  must  remain  in  fixed  relative  position,  the  reference  points  of  its  parts 
will  coincide  throughout  its  motion.  When  the  reference  point  is  in  C-space  obstacle 
X/Y,  the  following  constraint  is  in  eflfect:  “If  X  is  moving  and  Y  is  stationary,  a 
collision  has  occurred.”  In  other  words,  if  a  path  passes  through  X/F,  then  either 
X  must  be  stationary  or  Y  must  be  moving. 

Since  all  of  the  pairwise  C-spaces  were  computed  with  the  same  reference  point, 
we  now  superimpose  all  the  C-space  obstacles.  That  is,  all  the  obstacles  are  embed¬ 
ded  in  the  same  coordinate  system.  A  point  in  this  space  represents  a  displacement 
of  any  subset  of  the  parts  in  the  assembly,  and  the  C-space  obstacles  containing 
that  point  give  the  constraints  on  which  parts  may  be  members  of  the  moving  sub- 
assembly.  The  boundaries  of  the  obstacles  then  divide  the  plane  into  cells  such  that 
all  points  in  any  one  cell  are  in  the  same  set  of  C-space  obstacles.  The  set  of  super¬ 
imposed  C-spaces  is  the  interference  diagram  for  the  assembly.  Figure  6  shows  the 
interference  diagram  for  the  assembly  of  figure  5(a),  obtained  by  superimposing  all 
6  pairwise  C-spaces. 

Each  cell  in  the  interference  diagram  is  labeled  with  the  C-space  obstacles  that 
include  the  cell.  Each  label  is  a  constraint  on  any  subassembly  placed  in  that  cell. 
For  example,  the  cell  to  the  left  of  the  origin  in  figure  6  has  the  labels:  A/B^  A/C 
and  C/B,  If  the  reference  point  attached  to  A  is  placed  anywhere  in  that  cell,  A 
will  collide  with  B  and  C.  Similarly,  placing  C  in  that  cell  will  cause  a  collision  with 
part  B.  Importantly,  since  the  labels  B/A  and  B/C  are  missing  from  the  cell,  B 
can  safely  be  placed  there,  as  shown  in  figure  7(a).  Finally,  the  subassembly  {B^C} 
is  safe,  since  neither  B  nor  C  collides  with  the  remaining  objects  (figure  7(b)). 

The  constraints  in  a  cell  can  be  represented  and  analyzed  using  a  blocking 


11 


A/B 


Figure  6:  The  interference  diagram  for  the  assembly  of  figure  5(a).  The  small  labels 
indicate  the  C-space  obstacles  that  share  that  cell. 


(a)  (b) 

Figure  7 :  Two  subassemblies  that  may  be  placed  in  a  cell  of  the  interference  diagram. 
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(a)  (b) 


Figure  8:  Two  blocking  graphs  from  the  interference  diagram  of  figure  6. 

graph  [28].  A  blocking  graph  is  a  directed  graph  with  a  node  for  each  part;  a 
label  X/Y  in  the  cell  induces  a  directed  arc  from  node  X  to  node  Y  in  the  graph. 
If  label  X/Y  is  present,  part  Y  is  said  to  block  part  A'  in  that  cell.  If  a  proper 
subgraph  S  of  the  blocking  graph  has  no  outgoing  arcs  to  the  rest  of  the  graph, 
then  S  represents  a  subassembly  that  may  be  placed  in  the  cell,  and  S  is  called  free. 
Such  a  subgraph  exists  if  and  only  if  the  blocking  graph  is  not  strongly  connected.^ 
If  a  blocking  graph  is  not  strongly  connected,  then  one  of  its  strong  components  has 
no  outgoing  arcs  and  represents  a  subassembly  that  may  be  placed  in  the  cell.  The 
reduced  graph  of  strong  components  can  be  analyzed  in  a  straightforward  way  to 
identify  all  free  subassemblies  in  the  blocking  graph. 

Figure  8(a)  shows  the  blocking  graph  for  the  cell  we  considered  above.  It  is 
not  strongly  connected,  and  since  the  subassemblies  {B}  and  both  have  no 

outgoing  arcs,  each  may  be  placed  anywhere  in  the  cell  without  collision  (figure  7). 
Figure  8(b)  shows  the  blocking  graph  for  the  cell  just  above  the  origin  in  figure  6; 
it  is  strongly  connected,  so  no  proper  subassembly  can  be  placed  in  that  position. 

4.3  Removing  a  Subassembly 

When  a  subassembly  moves,  the  common  reference  point  moves  through  a  sequence 
of  cells  in  the  interference  diagram,  each  a  neighbor  of  the  previous  one.  Such  a 
sequence  of  connected  cells  is  called  a  path,  and  a  blocking  graph  can  be  associated 
with  it.  A  subassembly  that  can  follow  a  path  must  be  collision- free  at  each  point 
along  the  path,  so  the  blocking  graph  for  a  path  is  the  union  of  the  blocking  graphs 
associated  with  its  cells.  If  the  blocking  graph  for  the  path  is  not  strongly  connected, 
then  some  subassembly  can  follow  it  and  the  path  is  called  feasible. 

strongly  connected  component  {or  strong  component)  of  a  directed  graph  is  a  maximal  subset 
of  nodes  such  that  for  any  pair  of  nodes  (ni,n2)  in  this  subset,  a  path  connects  ni  to  712-  A  graph 
is  strongly  connected  if  it  has  only  one  strong  component. 
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Figure  9:  A  disassembly  path  for  the  assembly  in  figure  5(a)  and  the  blocking  graph 
it  induces.  Only  the  cells  of  the  interference  diagram  traversed  by  the  path  are 
labeled. 

Let  the  initial  cell  be  the  cell  containing  the  origin  of  the  interference  diagram. 
Similarly,  the  final  cell  is  the  free,  outermost  cell  of  the  interference  diagram;  any 
subassembly  in  this  cell  is  completely  separated  from  its  complement.  To  partition 
the  assembly,  a  feasible  path  must  be  found  that  connects  the  initial  cell  to  the  final 
cell.  Such  a  feasible  path  is  called  a  disassembly  path,  and  the  free  subassemblies  of 
the  corresponding  blocking  graph  may  be  removed  along  the  disassembly  path. 

Figure  9  shows  a  disassembly  path  for  the  assembly  of  figure  5(a),  and  its  blocking 
graph.  The  removable  subassembly  is  {B,C},  which  as  required  has  no  outgoing 
arcs  in  the  blocking  graph. 

The  search  for  a  disassembly  path  proceeds  by  extending  feasible  paths  from 
the  initial  region  outward.  The  initial  region  has  no  constraints,  since  in  any  valid 
assembly  all  parts  are  disjoint  in  their  initial  positions.  When  a  path  is  extended 
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from  a  region  to  one  of  its  neighbors,  the  constraints  for  any  C-space  obstacles  whose 
boundaries  are  crossed  are  added  to  the  blocking  graph  for  the  path.  If  at  any  point 
the  blocking  graph  becomes  strongly  connected,  the  path  is  discarded.  A  feasible 
path  that  reaches  the  final  region  is  a  disassembly  path,  and  the  free  subassemblies 
of  its  blocking  graph  can  be  removed  along  that  path. 

The  complexity  of  the  interference  diagram  itself  is  clearly  polynomial  in  the 
number  of  edges  of  the  parts,  and  the  diagram  can  be  constructed  in  similar  time. 
The  problem  is  that  the  number  of  paths  to  consider  might  be  very  large.  Not  aU  of 
these  paths  need  be  considered;  for  example,  when  two  paths  pass  through  the  same 
cell,  and  the  blocking  graph  of  one  path  is  a  subset  of  the  blocking  graph  of  the 
other,  the  more-constrained  path  may  be  discarded.  However,  such  simplifications 
do  not  reduce  the  complexity  of  the  problem,  which  is  exponential  in  the  general 
case. 

4.4  Generalizations 

The  interference  diagram  extends  directly  to  more  general  C-spaces  (hence,  mo¬ 
tions  of  the  subassemblies).  In  particular,  the  generalization  to  three-dimensional 
translations  of  polyhedra  is  clear;  the  interference  diagrain  is  an  xyz  C-space,  and 
the  regions  are  polyhedral.  This  is  quite  manageable  in  practice.  In  theory,  one 
can  carry  out  this  construction  in  arbitrary  C-spaces,  including  the  6-dimensional 
one  of  general  rigid  motion,  and  for  parts  with  arbitrary  curved  surfaces.  For  rigid 
motions,  the  reference  point  attached  to  each  part  is  replaced  by  a  reference  frame 
that  coincides  with  the  origin  in  the  initial  configuration.  A  motion  of  this  reference 
frame  can  then  refer  to  any  subset  of  the  parts.  As  before,  a  disassembly  path  is  a 
path  connecting  the  initial  ceU  (containing  the  original  placement  of  the  frame)  to 
the  final  cell  surrounding  the  rest  of  the  cells.  In  each  case  the  interference  diagram 
is  of  polynomial  complexity  and  can  be  constructed  in  polynomial  time.  However, 
implementing  such  constructions  is  difficult. 

To  allow  contacts  between  parts,  the  boundaries  of  the  C-space  obstacles  are 
considered  separate  regions,  in  which  the  tw'o  parts  do  not  collide.  Then  the  con¬ 
straint  X/Y  is  only  added  to  the  blocking  graph  for  a  path  as  that  path  is  extended 
into  the  interior  of  the  C-space  obstacle  XjY . 

While  exponential  in  the  worst  case,  a  partitioning  algorithm  for  arbitrary  disas¬ 
sembly  paths  that  searches  an  explicit  interference  diagram  may  well  prove  practical 
for  typical  industrial  assemblies.  We  have  not  implemented  the  method  to  test  this 
hypothesis;  instead,  the  next  sections  describe  special  cases  of  the  interference  dia¬ 
gram  that  yield  polynomial-time  partitioning  for  restricted  assembly  motions.  These 
cases  prove  very  useful  in  practice. 
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5  Partitioning  with  Single  Translations 

Very  few  assembly  operations  in  industry  require  the  kind  of  convoluted  paths  that 
the  general  partitioning  problem  allows;  in  fact,  an  assembly  that  needs  such  a 
complicated  motion  is  a  prime  candidate  for  redesign  [5].  Instead,  the  vast  major¬ 
ity  of  operations  can  be  accomplished  by  relatively  simple  motions,  such  as  single 
translations  to  infinity  or  extended  twisting  motions.  In  this  section  we  examine 
the  effect  on  the  partitioning  problem  w'hen  disassembly  paths  are  limited  to  single 
translations. 

A  single-translation  disassembly  path  is  a  ray  with  its  base  at  the  origin  of 
the  interference  diagram.  As  in  the  general  case,  the  cells  intersected  by  the  ray 
determine  a  set  of  constraints  on  subassemblies  that  may  be  removed  along  the  ray. 
A  ray  is  a  disassembly  path  if  the  blocking  graph  induced  by  these  constraints  is 
not  strongly  connected. 

Now  consider  what  happens  when  the  ray  rotates  to  new  disassembly  directions. 
The  blocking  graph  will  not  change  as  long  as  the  rotating  ray  intersects  the  same 
set  of  cells.  Indeed,  the  blocking  graph  changes  only  when  the  ray  begins  to  inter¬ 
sect  a  new  C-space  obstacle,  or  stops  intersecting  one  of  the  obstacles.  We  could 
partition  the  assembly  by  constructing  the  interference  diagram,  then  sweeping  a 
ray  through  it,  checking  the  blocking  graph  induced  at  each  change.  However,  the 
same  computation  can  be  performed  more  efficiently  by  considering  the  orientations 
of  the  ray  directly. 

For  planar  assemblies,  we  represent  the  set  of  all  disassembly  directions  by  the 
points  on  the  unit  circle  S^.  Projecting  each  C-space  obstacle  X/Y  onto  the  circle 
by  a  central  projection,  we  obtain  a  single  arc;  this  arc  consists  of  exactly  those 
directions  in  which  the  ray  wiU  intersect  X/Y?  For  example,  figure  10  shows  two 
polygonal  parts  from  figure  5(a),  their  C-space  obstacle,  and  the  projection  onto  the 
unit  circle.  The  C-space  obstacles  give  rise  to  n{n  —  1)  overlapping  arcs  on  5^;  the 
blocking  graph  only  changes  at  the  endpoints  of  these  arcs.  Checking  each  graph 
for  strong  connectedness,  a  disassembly  path  can  be  detected  if  one  exists. 

In  three  dimensions,  we  represent  the  disassembly  directions  as  points  on  the  unit 
sphere  S^.  Again  we  do  not  need  to  construct  the  interference  diagram  explicitly, 
but  only  its  projection  on  S^.  The  C-space  obstacle  for  two  polyhedral  parts  is  a 
polyhedron,  and  when  projected  onto  the  sphere  this  becomes  a  region  bounded  by 
arcs  of  great  circles  (see  figure  11).  This  region  is  the  set  of  all  translation  directions 
in  which  the  two  parts  will  collide.  The  edges  of  the  regions  for  aU  pairs  of  parts 
determine  an  arrangement  of  faces  on  5^;  the  blocking  graph  for  a  translational 

X  and  Y  aie  in  contact  in  their  initial  positions,  this  projection  must  be  handled  more 
carefully,  but  the  result  is  still  a  single  arc. 
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(a) 


(b) 


Figure  10:  (a)  Two  polygonal  parts  with  common  reference  point,  and  (b)  their 
C"Space  obstacle  projected  onto  the  unit  circle  (the  bold  arc). 


Figure  11:  (a)  Two  simple  polyhedral  parts,  and  (b)  their  C-space  obstacle,  pro¬ 
jected  onto  the  sphere. 


disassembly  path  only  changes  at  the  boundaries  of  these  regions.  Each  face  is 
labeled  with  the  labels  that  project  onto  it.  In  practice,  the  C-space  obstacles  need 
not  be  constructed  explicitly;  their  projections  can  be  built  by  triangulating  the 
faces  of  the  polyhedra,  computing  and  projecting  C-space  obstacles  for  each  pair  of 
triangles,  and  taking  their  union  on  the  sphere. 

For  an  assembly  of  n  polyhedra  with  a  total  of  v  vertices,  the  n{n  -  1)  regions 
on  5^  are  bounded  by  a  total  of  arcs.  The  arrangement  they  determine  has 

complexity  and  can  be  calculated  in  the  same  time  bound  [6].  The  blocking 

graphs  can  be  found  by  traversing  the  arrangement,  making  small  changes  as  each 
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cell  boundary  is  crossed.  Checking  a  blocking  graph  for  strong  connectedness  may 
require  0{v?)  time,  for  a  total  of  0{n^v^)  time  to  partition  a  polyhedral  assembly 
with  single  translations.  If  all  partitionings  are  desired,  they  can  be  found  by  ana¬ 
lyzing  the  reduced  graph  of  strong  components  for  each  blocking  graph.  To  find  an 
assembly  sequence  of  operations  with  single  translations  requires  n  —  1  applications 
of  the  partitioning  algorithm,  or  0{n^v'^)  time. 

6  Partitioning  with  Infinitesimal  Rigid  Motions 

This  section  describes  another  special  case  of  the  interference  diagram  that  is  useful 
when  an  assembly  cannot  be  partitioned  using  single  translations,  or  when  other 
partitionings  are  desired  to  optimize  the  assembly  sequence  in  some  way.  It  is  moti¬ 
vated  by  the  following  observation:  any  disassembly  path,  no  matter  how  convoluted, 
must  begin  with  the  initial  cell  (the  origin  of  the  interference  diagram)  followed  by 
one  of  its  neighbors.  Hence  a  necessary  condition  on  removable  subassemblies  can 
be  computed  by  exploring  the  local  topology  of  the  neighborhood  of  the  initial  cell 
in  the  interference  diagram.  The  subassemblies  that  can  follow  a  path  from  the 
initial  cell  to  one  of  its  immediate  neighbors  are  exactly  those  that  can  move  a  very 
small  distance  in  the  assembly.  Each  subassembly  that  satisfies  this  local  freedom 
constraint  can  then  be  passed  to  a  path  planner  or  other  evaluation  technique. 

The  topology  of  the  initial  cell’s  neighborhood  is  determined  wholly  by  the  con¬ 
tacts  between  parts  in  their  initial  placement;  no  other  constraints  come  into  play 
until  a  subassembly  moves  a  finite  distance.  Two  parts  are  in  contact  when  their 
boundaries  intersect.  In  experiments,  real-world  assemblies  have  proven  to  be  quite 
constrained  by  their  contacts  [3,  16,  30],  resulting  in  very  few  subassemblies  that 
must  be  further  tested  for  removability. 

Let  us  consider  disassembly  motions  consisting  of  arbitrary  rigid  motions  (com¬ 
bining  translation  and  rotation)  in  3D  for  polyhedral  parts;  the  2D  and  translation 
cases  are  just  simpler.  The  interference  diagram  in  this  case  is  a  6D  C-space  con¬ 
struction.  The  neighborhood  of  interest  is  an  infinitesimal  5D  sphere  surrounding 
the  origin  of  the  interference  diagram.  This  sphere  is  cut  into  faces  of  dimension 
0  (a  vertex)  to  5  by  the  surfaces  of  the  interference  diagram  that  neighbor  the  ini¬ 
tial  cell,  and  each  face  is  labeled  with  the  blocking  graph  of  the  cell  containing  it. 
By  checking  these  blocking  graphs  for  strong  connectedness,  we  can  determine  the 
subassemblies  that  are  free  to  move  a  small  distance  in  the  assembly. 

We  represent  the  infinitesimal  sphere  by  the  unit  sphere  5^.  We  could  construct 
this  arrangement  on  the  sphere  by  building  the  6D  interference  diagram  for  rigid 
motions,  then  cutting  5^  by  the  tangent  hyperplanes  to  surfaces  through  the  origin. 
However,  the  tangent  hyperplanes  and  arrangement  on  can  be  computed  more 
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Figure  12:  A  point-plane  contact. 


directly  from  the  contacts  of  the  assembly. 

The  constraints  that  contacts  impose  on  the  relative  motion  of  the  two  contacting 
bodies  have  been  studied  in  depth  (see  for  instance  [22,  24,  26]).  An  infinitesimal 
rigid  motion  can  be  described  by  a  6D  vector  dX  =  (dx,  dy.  dz,  da,  d/3,  d'/)  giving  the 
components  of  translation  and  rates  of  infinitesimal  rotation  about  the  coordinate 
axes,  respectively.  Consider  two  parts  P,  and  Pj  in  contact,  and  let  V  be  a  vertex 
of  Pi  (figure  12).  A  motion  dX  applied  to  P,  causes  V  to  undergo  a  translation 
t  =  JvdX,  where  Jv  is  a  constant  3x6  Jacobian  matrix  relating  the  motion  of 
Pi  to  the  translation  of  V.  Assume  that  P,  and  Pj  are  in  contact  such  that  V  is 
contained  in  a  face  of  Pj  with  outward  normal  n.  Then  dX  causes  V  to  penetrate 
Pj  when  n  •  t  <  0,  to  break  the  contact  when  n  -t  >  0,  and  to  slide  in  F  when 
n-t  =  0.  The  hyperplane  n  •  t  =  0  (the  sliding  case)  is  the  tangent  hyperplane  at 
the  origin  to  a  surface  of  the  interference  diagram.  Hence  the  motions  dX  allowed 
by  the  contact  are  those  satisfying  n  •  (JydX)  >  0. 

The  motions  allowed  by  other  contacts  between  polyhedra  can  be  expressed  with 
conjunctions  of  point-plane  contact  constraints.  For  instance,  consider  a  contact 
between  a  face  of  P,  and  a  face  of  Pj  with  outward  normal  n  (see  figure  13).  The  set 
of  motions  dX  allowed  by  this  contact  is  the  intersection  of  the  closed  half-spaces 
n  •  (Jv^dX)  >  0  computed  for  each  of  the  vertices  14  of  the  convex  hull  of  the 
intersection  of  the  two  faces.  In  figure  13  the  vertices  14  are  circled.  See  [13,  28]  for 
the  constraints  arising  from  other  polyhedral  contact  types. 

The  intersection  of  the  half-spaces  of  point-plane  constraints  between  each  pair 
of  parts  Pi  and  Pj  is  a  closed  convex  6D  polytope  in  the  space  of  infinitesimal 
motions.  Pj  blocks  the  motion  of  P,-  for  all  motions  contained  in  this  polytope,  and 
for  no  motions  outside  the  polytope.  The  intersection  of  these  polytopes  with  the 
sphere  5®  determines  an  arrangement  of  regions  of  dimensions  0, . . . ,  5  on  5^.  These 
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Figure  13:  Plane-plane  contact  expressed  as  point-plane  constraints. 

regions  capture  the  local  topology  of  the  origin  of  the  interference  diagram. 

Let  there  be  m  <  (”)  pairs  of  parts  in  contact  in  the  assembly,  described  by 
a  total  of  c  €  O(v^)  point-plane  constraints,  where  c  >  m.^  Then  the  arrange¬ 
ment  of  poly  topes  on  5®  has  O(m^c^)  cells  [2],  and  can  be  constructed  by  finding 
the  arrangement  of  hyperplanes  in  O(c^)  time  using  a  topological  sweep  [8]  then 
merging  cells.  Since  each  cell’s  blocking  graph  can  have  at  most  2m  edges,  checking 
it  for  strong  connectedness  can  be  performed  in  time  0{m).  Hence  identifying  a 
subassembly  that  is  locally  free  has  total  running  time  of  0{c^  4-  m‘*c^).  This  bound 
wiU  improve  if  and  when  better  algorithms  to  calculate  arrangements  in  high  di¬ 
mensions  are  developed.  Note  that  in  practice,  many  free  subassemblies  wiU  need  to 
be  generated  from  the  reduced  graph  of  strong  components  for  each  blocking  graph, 
since  some  locally-free  subassemblies  may  not  satisfy  global  constraints. 

A  version  of  this  algorithm  has  been  implemented  and  tested  on  a  number  of  in¬ 
dustrial  assemblies.  The  experimental  system  allows  parts  with  planar,  cylindrical, 
and  some  helicoidal  surfaces.  The  implemented  version  of  the  algorithm  considers 
the  set  of  all  infinitesimal  translations  dX  =  {dx,dy,dz),  plus  a  finite  set  of  “sug¬ 
gested  motions”  inferred  from  nonplanar  contacts  in  the  assembly.  For  example, 
a  cylindrical  contact  suggests  pure  rotation  about  its  axis;  a  threaded  contact  be¬ 
tween  two  helicoidal  surfaces  suggests  a  screwing  motion;  etc.  The  set  of  suggested 
generalized  motions  is  incomplete  but  accounts  for  most  motions  required  by  ac¬ 
tual  assemblies.  Figure  14  shows  a  model-aircraft  combustion  engine  with  42  parts; 
the  planner  identifies  all  33  subassemblies  of  the  engine  that  satisfy  local  freedom 

®In  most  real  assemblies,  m  ss  n,  c  and  c  m. 
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Figure  14:  Model- aircraft  engine 


constraints  (of  which  25  satisfy  global  constraints)  in  8.5  seconds  on  a  DECstation 
5000/200. 

7  Discussion 

We  have  shown  that  the  partitioning  problem  is  NP-complete  for  assemblies  of 
simple  polygons  in  the  plane  when  arbitrary  motions  are  allowed  to  separate  the 
two  subassemblies.  On  the  other  hand,  when  the  motions  allowed  for  assembly 
are  limited  to  single  translations  or  infinitesimal  rigid  motions,  partitioning  is  in  P. 
This  raises  the  question:  under  what  conditions  can  partitioning  be  accomplished 
in  polynomial  time? 

One  factor  appears  to  be  the  “complexity”  of  the  motions  allowed  to  separate 
the  two  subassemblies.  For  instance,  one  way  to  measure  the  complexity  of  a  trans¬ 
lational  path  is  by  the  number  of  straight-line  translations  it  includes.  In  separate 
work,  Halperin  and  Wilson  [12]  investigate  the  problem  of  partitioning  an  assembly 
with  a  small  number  t  of  translations.  The  case  f  =  1  is  the  single  translation  case 
described  above.  For  the  case  t  =  2,  a  disassembly  path  (in  two  dimensions)  can 
be  described  by  a  triplet  {x,y,<f>),  where  (x,y)  is  the  endpoint  of  the  first  transla¬ 
tion,  and  </>  is  the  orientation  of  the  second  translation  to  infinity.  For  instance,  the 
disassembly  path  shown  in  figure  9  is  of  this  type. 

As  in  the  single  translation  case,  the  cells  of  the  interference  diagram  cut  by  a 
two-translation  path  (a;,  y,4>)  do  not  change  when  x,  t/,  and  <f>  vary  slightly.  The  set  of 
cells  cut  by  the  path  changes  only  at  certain  critical  values.  These  values  determine 
surfaces  in  (x,j/,^)-space,  dividing  this  space  into  regions;  each  region  represents 
a  set  of  two-translation  paths  that  have  the  same  blocking  graph.  Checking  the 


graphs  for  strong  connectedness  yields  a  two-translation  partitioning  algorithm. 

In  theory,  this  easily  extends  to  disassembly  paths  of  more  than  two  translations. 
A  t-translation  path  is  described  by  a  tuple  (xi,  j/i, . . . where  (xj,  j/,) 
is  the  endpoint  of  the  i-th  translation  in  the  path.  The  cells  of  the  interference 
diagram  now  determine  hyper  surfaces  in  a  (2t—l)- dimensional  space  of  disassembly 
paths.  For  any  fixed  t  this  approach  yields  a  polynomial  time  partitioning  algorithm; 
however,  the  number  of  regions  in  the  arrangement  can  be  exponential  in  t,  so  the 
resulting  algorithm  will  be  practical  only  for  very  small  values  of  t. 

In  experimental  assembly  planning  systems,  an  additional  constraint  is  often 
added  to  the  partitioning  problem,  requiring  that  both  5  and  A  \  5  be  connected. 
As  mentioned  at  the  end  of  section  3,  a  subassembly  is  considered  connected  if  the 
union  of  its  parts  is  a  connected  set.  This  constraint  is  useful  in  practice,  since  un¬ 
connected  subassemblies  are  usually  more  difficult  to  grasp,  fixture,  and  manipulate. 
When  the  connectedness  constraint  is  included,  partitioning  assemblies  of  polyhe- 
dra  is  NP-complete.  However,  it  remains  an  open  problem  whether  partitioning 
an  assembly  of  simple  polygons  into  two  connected  subassemblies  is  NP-complete. 
While  partitioning  in  2D  is  less  applicable  to  real-world  assemblies  than  the  3D  case, 
solving  this  open  problem  might  lend  insight  into  the  assembly  sequencing  problem. 

Another  open  problem  regards  checking  many  blocking  graphs  for  strong  con¬ 
nectedness.  When  disassembly  motions  are  limited  to  single  translations  or  in¬ 
finitesimal  motions,  the  interference  diagram  becomes  an  arrangement  of  regions, 
each  associated  with  a  blocking  graph.  The  blocking  graphs  of  neighbor  regions  dif¬ 
fer  only  slightly:  if  a  single  surface  separates  two  regions,  then  the  blocking  graphs 
for  the  regions  differ  by  at  most  one  arc.  Checking  whether  a  directed  graph  is 
strongly  connected  requires  time  linear  in  the  number  of  edges,  so  the  time  to  check 
all  the  blocking  graphs  dominates  the  running  time  of  each  algorithm.  Reducing  the 
cost  of  these  closely-related  connectedness  checks  would  reduce  the  running  time  of 
the  algorithms  significantly.  However,  at  this  time  no  improvement  over  the  naive 
method  has  been  found. 

Finally,  the  techniques  considered  in  this  paper  only  consider  the  constraints 
arising  directly  from  the  geometry  of  the  parts.  There  are  many  other  constraints 
on  assembly  plans,  such  as  the  need  to  grasp,  fixture,  and  manipulate  the  assemblies. 
Constraints  also  arise  from  the  effects  of  uncertainty  or  flexible  parts  on  the  assembly 
process.  Much  work  remains  to  be  done  in  these  areas. 


Appendix 

We  discuss  here  the  choice  of  x,’s  that  are  used  in  the  assignment  mechanism  of  our 
construction  of  Section  3.  The  Xj’s  indicate  the  relative  horizontal  positions  of  the 
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parts  of  the  moving  subassembly  S  and  are  preserved  throughout  the  motion  of  S. 

Assume  that  S  =  where  Rq  denotes  the  key  and  Ri  denotes 

either  Ui  OTUi,l<i<  n.  We  say  that  Ri  is  at  position  x  and  write  pos{Ri)  =  x  iff 
the  x-coordinate  of  the  bottom-left  vertex  of  the  rectangle  Ri  is  x. 

Lemma:  Let  Xi  —  10  •  a,-,  i  =  0,1, . .  .,n,  where  aQ,ai,  —  ,a„  is  an  integer  se¬ 

quence  such  that  all  pairwise  differences  a,-  -  aj  are  distinct  when  i  ^  j.  When 
|xp  -  pos(ili)|  <  4,  p  7^  i,  then  for  all  Rj,  j  ^  i,  we  have  that  |x,  -  po${Rj)\  >  4. 

Proof  At  some  time  during  the  motion  of  5,  Ri  is  4-close  to  the  position  Xp,  that 
is  |xp  -  pos(il;)|  <  4.  Consider  now  any  Rj,  j  ^  i,  and  any  posistion  x,  with  q  7^  j. 
Let  us  bound  the  quantity  \xq-pos{Rj)\  from  below  using  the  triangular  inequality: 

jx,  -  pos{Rj)\  >  |(xp  -  pos{Ri))  -  {xq  -  pos{Rj))\  -  |xp  - 

or  equivalently, 

|xg  -  pos(J?j)|  >  |(Xp  -  Xq)  +  {pos{Rj)  -  pos(iEi))|  -  |xp  -  pos(Ji:i)| 

But  since  |xp  —  po5(iZ,)l  <  4  and  pos{Rj)  —  pos{Ri)  =  Xj  —  x,-  =  10  •  aj  —  10  •  c,'  the 
above  inequality  becomes: 

|xg  -  pos{Rj)\  >  10  •  jcp  -  ttq-  {aj  -  a,)|  -  4. 

From  the  hypothesis  about  an  the  quantity  |ap-ag-l-(aj-a,)|  =  |(aj-ag)-(aj-ap)| 
must  be  at  least  one,  since  j  ^  i  and  i  ^  p.  Thus  we  get: 

\pos{Rj)  —  Xgl  >  10  •  1  —  4  >  4. 

which  completes  the  proof.  O 

The  above  Lemma  guarantees  that  for  any  k,  w'hen  Rk  is  close  to  position  Xp, 
p  k,  and  needs  to  go  through  the  hole  that  has  been  created  at  this  position, 
none  of  the  Rq,Ri,...,  ...Rn  is  close  to  a  position  x,.  That  is,  all  of 

Ro,  Ri, . . . ,  Rk-i,  Rk+i,  •  •  ■  Rn  are  in  the  wide  parts  of  the  assignment  mechanism 
and  can  follow  the  constrained  motion  of  Rk- 

The  key  element  in  the  above  proof  is  the  property  of  the  sequence  a,-,  namely 
that  the  pairwise  diffrences  of  its  terms  are  all  distinct.  The  10  is  just  a  scaling 
factor  that  gives  extra  space  for  the  width  of  the  parts.  It  is  not  hard  to  choose  the 
a,’s;  a  straightforward  example  of  sequence  that  has  all  its  terms  distinct  is  ai  =  2*. 
How'ever,  this  sequence  yields  an  exponential  length  for  the  assignment  construct. 
Using  a  result  of  Erdos  [11],  it  is  possible  to  select  the  Oi’s  in  such  a  way  that 
max{ao, . .  .,a„}  =  O(n^).  Then  the  physical  length  of  the  assignment  mechanism 
becomes  a  polynomial  of  n. 
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